Swarming ring patterns in bacterial colonies exposed to ultraviolet radiation 
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We report a novel morphological transition in a Bacillus subtilis colony initially growing under 
ambient conditions, after ultraviolet radiation exposure. The bacteria in the central regions of the 
colonies are observed to migrate towards the colony edge forming a ring during uniform spatial 
exposure. When the radiation is switched off, the colonies were observed to grow both inward into 
the evacuated regions as well as outward indicating that the pattern is not formed due to depletion 
of nutrients at the center of the colony. We also propose a reaction-diffusion model in which waste- 
limited chemotaxis initiated by the UV radiation leads to the observed phenomenology. 

PACS numbers: 87.18.Hf, 47.20.Hw, 87.50. Gi, 87.23. Cc 



Bacterial colony growth on semi-solid agar has served 
as a model system to understand self-organization of 
microscopic biological elements to form complex pat- 
terns [Q-pl. Although interaction between individuals 
may be well known, it is difficult to determine how com- 
plex patterns form for example due to chemotaxis. Ran- 
dom walk and reaction diffusion models have been also 
developed to find the minimal mechanisms that lead to 
the formation of the observed patterns J3|-f7j . In most of 
these studies, changes in the environment occur due to 
the slow depletion of nutrients and excretion of waste. 

The modeling of the response of living organisms to a 
change in environment is an important issue of current 
interest. An example is the effect of ultraviolet (UV) ra- 
diation on biological systems prompted by stratospheric 
ozone depletion. Recent studies have modeled the re- 
sponse to externally imposed change on bacterial colonies 
using the Fisher wave equation with added convection 
terms [^|J^]. Experiments using ultraviolet illumination 
pattern changes on Bacillus subtilis colonies have shown 
the relevance of such modeling [ ^0| . However, a delay in 
the response of the bacterial colony to changes in the en- 
vironment was also noted. These studies raise interesting 
questions on how to model biological systems that may 
undergo growth mode changes and show history depen- 
dence in response to environmental changes. 

In this paper, we report the spatio-temporal response 
of colonies of Bacillus Subtilis growing on nutrient rich 
agar to a sudden change in the environment. Bacillus 
Subtilis is a 9/im x 2/xm motile bacterium that can be 
often found in soil. The change is hostile and is imple- 
mented by exposing the colony to uniform UV radiation. 
The colony which grows uniformly before exposure to 
UV radiation, is not observed to simply stop due to the 
radiation but undergo a morphology change. Most of 
the bacteria are observed to migrate to the edge of the 
colony even though this does not help the bacteria to 
evade the radiation. To systematically study the spatio- 



temporal response of the colony, we measured the change 
in the bacterial density as a function of the strength of the 
UV radiation for the wild type and DNA repair mutants, 
uvrA and uvrA recA of the bacterium. We also propose a 
simple reaction-diffusion model which reproduces most of 
the observed phenomenology. The model operates with 
three densities: bacteria, waste, and chemoattractant, 
and describes a waste-limited chemotaxis initiated by the 
UV radiation. The bacteria migrate towards the edges 
because the waste concentration is reduced at the edges 
of the colony. 

Experiments were performed in 15 cm diameter Plex- 
iglas petri dishes containing a thin layer of nutrient agar 
(7 grams/liter of bacto peptone and 3 grams/liter agar). 
Under these conditions the growth of the bacterial colony 
can be approximately modeled by the Fisher wave equa- 
tion jll],[l0). The bacterial colony was initiated with an 
inoculating needle at a single point source at the center of 
the petri dish and the bacteria grew in a thin layer at the 
surface of the medium. The colony is first allowed to grow 
in the absence of UV radiation. The initial adaptation 
time and subsequent growth rate of the colony depends 
on the temperature and for most of the experiments a 
temperature of 24°C was used. At this temperature, the 
colony front was observed to grow slowly for the first 26 
hours and then increased reaching a constant front speed 
Vf = 0.3/im/s. The colony was then exposed to uniform 
UV radiation over a time interval t uv using a fluores- 
cent bulb, (Philips F15T8) which radiates in the range 
of 310-400 nm (UV-A and UV-B). An example of the 
spatio-temporal response in a colony of Bacillus subtilis, 
uvrA under UV radiation is shown in Fig. [l](a)-(c). 

We note that the circular symmetry of the patterns is a 
consequence of the point inoculation and the subsequent 
disk shaped colony which develops before the colony is 
exposed to UV radiation. If the bacteria is grown in a 
ribbon shape using a line inoculation pattern, they still 
migrate to the edges. For simplicity, we performed quan- 
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titative experiments with colonies using point inocula- 
tion. 

To further investigate the cause of the migration of 
the bacterial population, we observed the recovery of the 
colony after exposure to radiation. Images of the recover- 
ing colony are shown in Fig. |l](d)-(f) and it can be noted 
that the colony recovers almost uniformly in both direc- 
tions. The colony grows outward to fresh regions and 
also inward into the evacuated central regions. 

The radius of the circular colony before exposure, dur- 
ing exposure and after exposure was plotted as a func- 
tion of time in Fig. ^. We observed that the growth of 
the front of the colony during exposure is negligible com- 
pared to the growth before exposure. After radiation is 
switched off, it can be observed that cell division resumes 
following a "recovery" time t r , which increases with the 
strength of the radiation. Although, the front speed on 
the outside edge is greater than the inside edge, clearly 
the central regions of the colony are hospitable. There- 
fore, the initial evacuation of the colony at the center is 
not due to depletion of nutrients. Thus the formation 
of the swarm ring like patterns has clearly a different 
mechanism than observed previously 0,0- 

To quantitatively show the rearrangement of the bac- 
terial colony as a function of time, the bacterial density 
n(r) as a function of distance r from the center is plotted 
in Fig. ||(a). The density of the bacteria scales linearly 
with light intensity and therefore the density distribution 
of the colony is obtained from the images using a CCD 
camera. At approximately 24 hours the density of the 
bacteria is considerably higher at the edge of the colony 
compared to the colony center. After 70 hours the bacte- 
ria almost completely evacuate the center of the colony. 

In Fig. ||(b), we show the dependence of the strength 
of the UV radiation / on the colony. The cross section of 
the colony was plotted after the colonies were allowed to 
grow and were exposed to UV radiation for the same time 
duration t uv . The front of the colony was also observed 
to advance faster with respect to the higher intensity, but 
was still much slower compared to the v f before exposure 
to UV light (see Table 1). The total population of the 
colony as a function of time was measured by integrat- 
ing the population densities and is shown in Fig. ||(c). 
For lower /, the population when exposed to radiation 
increases at a small rate initially, and then decreases. 
However, for higher / the total population always de- 
creases. 

Qualitatively similar phenomena was observed in the 
two other strains of Bacillus subtilis. Although the den- 
sity was observed to increase in the wild type strain near 
the edge of the colony, the center of the colony was not 
observed to vacate completely. At the same UV radi- 
ation strength, the total population and the front were 
also observed to grow at a faster rate compared to uvrA. 
In the case of uvrA recA, the colony was observed to 
have a higher density at the edges but the total popu- 



lation decreased strongly with time. Thus the behav- 
ior of the recA strain can be thought of as intermediate 
between that observed in the wild and the uvrA recA 
strains. This is consistent with the fact that while all 
mechanisms of DNA repair are present in the wild type, 
one mechanism of repair is absent in uvrA, and two are 
absent in the urvA recA strain Jl3|-^6|. Therefore the 
rate and strength of the density distribution appears to 
be related to the efficiency of the DNA repair mechanism. 

In an attempt to elucidate the mechanism of this mor- 
phological transition, we propose the following simple 



model |17|. The model consists of three reaction-diffusion 
equations for bacteria concentration 6(r,i), waste con- 
centration w(r,t), and chemoattractant concentration 
c(r,i): 

8 t b = fb(l -w)- n(l - f)b - V • J c + D b V 2 b (1) 
d t w = fb(l -w) + D w V 2 w (2) 
Btc = (1 - f)b(l - w) + D c V 2 c (3) 

The first terms in the r.h.s. of Eqs.(|T]),(||) describe 
the bacterial growth and the accompanying waste pro- 
duction, respectively. The bacterial growth is limited by 
the local concentration of waste w and saturates when 
waste concentration approaches 1. Constant / = 1 in 
the absence of the UV radiation and becomes small when 
UV radiation is switched on. According to Eq.(||), when 
/ < 1 , bacteria emit chemoattractant c at the rate which 
is proportional to the bacteria concentration and also is 
limited by the waste concentration. The term — — f)b 
in Eq.(|l|) is responsible for the slow destruction of bac- 
teria by the UV radiation observed experimentally. The 
chemotactic flux J c = ab(l — b)Vc is directed towards the 
gradient of the concentration of chemical c, and it satu- 
rates at large bacterial concentration b — ► 1 (hard core 
repulsion). The last terms in Eqs.(0)-(||) describe linear 
diffusion of the components. For simplicity, we assume 
that the diffusion constants of all three fields are equal 
and rescale them to unity. 

The initial colony expansion without UV (/ = 1) is de- 
scribed by the first two equations with J c = 0. Asymp- 
totically w — b — > 0, and the system is reduced to a single 
Kolmogorov-Petrovsky-Piskunov-Fisher equation ff§] 



d t b = /6(I-&) + V 2 &, 



(4) 



which has a propagating front solution with the front 
speed Vf = 2. 

The simplest way to model the effect of UV radiation 
is to set / = at some time t = to. The growth of the 
bacteria is arrested, and the colony expansion stops. The 
subsequent dynamics of the population is determined by 
the interplay of diffusion and chemotaxis. In the model, 
the bacteria immediately begin to emit chemoattractant 
c in the regions where waste w is below the limiting level 
w = 1, i.e. near the edges of the colony. A ring of c is 
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formed, and the exponentially decaying tails of this ring 
extend towards the center of the colony thereby pro- 
voking the bacteria to concentrate near the edges of the 
colony. 

We integrated the model (jl])-(||) numerically in a 
square 150x150 with periodic boundary conditions start- 
ing from a small radially-symmetric non-zero perturba- 
tion of b in the center for /i = 0.02, a — 2. Variations in 
the strength of the UV radiation can be simulated by us- 
ing different small values of f(t > to) instead of zero. In a 
typical simulation we switched / from 1 to 0.1 at t = 20, 
when the radius of the colony was r » 21. The resulting 
ring pattern and the radial density distribution as a func- 
tion of time are shown in Fig. ||. As in the experiment, 
after the UV radiation is turned on, the population starts 
to assemble near the edge (cf. Fig. ||(a)). Due to slow 
bacterial growth, the outer front continues to expand, 
however, at a much slower rate. The total population 
size slightly grows initially at t > to, but later declines. 
All this agrees completely with observed phenomenology. 
A more detailed discussion of the model and its compar- 
ison with experiments will be reported elsewhere. 

In summary, colonies of Bacillus subtilis are observed 
to show a swarming ring like pattern when exposed to 
UV radiation. Thus the mode in which the bacterial 
colonies grow after exposure to UV radiation is very dif- 
ferent from the Fisher wave like mode in which it grows 
prior to exposure. Using reaction-diffusion modeling, we 
have shown that the migration of colony under UV ra- 
diation is related to chemotaxis which is limited by the 
waste that Bacillus subtilis produces. 
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v g (UV) 
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30 


0.3 0.004 


0.016-0.020 


12 


0.3 0.006 


0.013-0.016 


7 


0.3 0.010 


0.011-0.013 



TABLE I. The front and group velocity of Bacillus subtilis 
uvrA as a function of intensity / of the UV radiation. 



Corresponding author 
[1] M. Matsushita and H. Fujikawa, Physica A 168, 498 
(1990). 

[2] E. O. Budrene and H. C. Berg, Nature 349, 630 (1991). 
[3] E. Ben- Jacob, et al, Nature 368, 46 (1994) . 
[4] L. Tsimring, et al, Phys. Rev. Lett. 75, 1859 (1995). 
[5] M. P. Brenner, L. S. Levitov, and E. O. Budrene, Bio- 

phys. J. 74, 1677 (1998). 
[6] Y. Kozlovsky, et al, Phys. Rev. E 59, 7025 (1999). 
[7] R. Tyson, S. R. Lubkin, and J. D. Murray, Proc. R. Soc. 

Lond. B 266, 299 (1999). 
[8] D. R. Nelson and N. M. Shnerb, Phys. Rev. E 58, 1383 

(1998). 



3 



FIG. 1. (a)-(c) The images of Bacillus subtilis, uvrA as a 
function of time t uv under ultraviolet (UV) radiation, (a) t uv 
= h (b) tuv = 25 h and (c) t uv = 45 h. The UV radiation 
of strength I — 30W/m 2 is switched on after the colony was 
allowed to grow for time t — 50 hours. The bacterial popu- 
lation is observed to migrate to the edge of the colony. The 
bright point at the center of the colony corresponds to the 
point of inoculation, (d)-(f) Recovery images of the bacterial 
population as a function of time t r after the UV radiation is 
switched off. The bacterial colony was allowed to recover for 
(d) tr = 30 h (e) t r = 40 h and (f) t r = 55 h. The bacterial 
population is observed to grow both outward as well as into 
the vacated inner region. 



Initial growth 

J 


UV 




ooccooa 


^^^^^ 



40 



80 
t (hour) 



120 



160 



FIG. 2. The radius of the colony as a function of time from 
the inoculation (7 = 30W/m 2 ). The outer edge of the colony 
is plotted using (•), and the inner edge using (o). The bacte- 
rial colony is observed to recover both inwards and outwards 
in the absence of UV radiation. The dashed line corresponds 
to the point where population density has decreased to less 
than 10 % of the maximum value at the edge. 
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FIG. 3. (a) The population density n{r)/n ma x as a func- 
tion of distance r from the inoculation point and as a func- 
tion of time for I = 12W/m 2 . n max is the maximum density 
at tuv = 0. (The data is shifted by approximately 2 mm 
from the inoculation point to avoid the puncture mark.) (b) 
The n(r)/n ma x as a function of radiation strength I at for 
46 hours. The ring is observed to form faster as the intensity 
of UV radiation is increased. The colony front for the three 
cases have been offset to align at the colony edge for clarity, 
(c) The fractional change of the total population as a func- 
tion of time after the UV radiation is turned on for the three 
intensities. The dashed line is the fractional increase in the 
population in the absence of UV radiation. 



FIG. 4. The bacterial density b(r) as a function of distance 
r from the center at 10 unit time intervals after the radiation 
was effectively turned on in Eqs.(|l|)-(|^) for /i = 0.02, a = 2. 
Inset: A ring formed using the proposed model is similar to 
that observed in the experiments. 
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